Benchmark calculations for ^H, ''He, ^^O and ^°Ca with ab-initio coupled-cluster theory 
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We present ab-initio calculations for ''H, ^He, ^'^O, and ""^Ca based on two-nucleon low-momentum 
interactions Viowfc within coupled-cluster theory. For and ^He, our results are within 70keV and 
lOkeV of the corresponding Faddeev and Faddeev-Yakubovsky energies. We study in detail the 
convergence with respect to the size of the model space and the single-particle basis. For the heavier 
nuclei, we report practically converged binding energies and compare with other approaches. 
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I. INTRODUCTION 

Ab-initio few- and many-body methods have been used 
with great success to explore the structure of light nu- 
clei based on microscopic two- and three-nucleon in- 
teractions. For nuclei with A < 12 nucleons, several 
techniques provide practically exact solutions to the nu- 
clear many-body problem and have been benchmarked to 
agree within numerical uncertainties for the ^He ground- 
state energy and radius obtained from the nucleon- 
nucleon (NN) Argonne potential These meth- 
ods include Faddeev-Yakubovsky equations 0], varia- 
tional approaches @,S,l5|,|6i, the Green's function Monte 
Carlo (GFMC) method Q, the No-Core SheU Model 
(NCSM) [si, ^ , and the effective interaction hyperspher- 
ical harmonics method [13, [HI • Many of the above ap- 
proaches are, however, restricted to the lightest nuclei. 

In recent years, the coupled-cluster method [l3,[i3| was 
reintroduced in nuclear physics as a tool for ab-initio nu- 
clear structure calculations [H, [H, [H [13, [3, HI, H^, HH . 
Coupled-cluster theory is size-extensive and scales rather 
gently with an increasing number of nucleons and with 
the size of the model space, and therefore has the poten- 
tial to extend the reach to medium-mass nuclei. In this 
paper, we address the question of whether the coupled- 
cluster method is as precise for few-body systems as the 
well-established methods. Several findings suggest that 
this is indeed the case. Mihaila and Heisenberg per- 
formed a microscopic calculation of the electron scatter- 
ing structure function for ^^O and found excellent agree- 
ment with experimental data. Their calculations are 
based on a particle-hole energy expansion of the cluster 
operator. More recent applications [H, [H, [llllS, IHI 
follow the "standard" approach of coupled-cluster calcu- 
lations from quantum chemistry [13, [3, [2^, [1^ . In these 
calculations, several results for helium isotopes [13, [131 
were found to be in good agreement with exact diago- 
nalizations in sufficiently small model spaces and with 
corresponding renormalized interactions. However, ex- 
cept in the recent study of three-nucleon forces {3NF) in 
coupled-cluster theory f21'| , this approach has never been 



compared in detail to well-established few-body methods 
in the larger model spaces that are needed for conver- 
gence with modern NN interactions. It is the purpose of 
this work to fill this gap in nuclear physics, and to place 
the coupled-cluster method in the group of ab-initio ap- 
proaches. 

This paper is organized as follows. In Section|lll we be- 
gin with a brief discussion of the coupled-cluster method, 
of the low-momentum interactions and the employed ba- 
sis spaces. Our main results for and ^He are presented 
in Section inn and for ^^O and ^^Ca in Section [iVl We 
conclude with a summary in Section |Vl 



II. METHOD, INTERACTIONS, AND MODEL 
SPACES 

A. Coupled-cluster method 

Coupled-cluster theory was invented by Coester and 
Kiimmel almost fifty years ago [13, [i3- During the 1970s, 
this approach was further developed and found many ap- 
plications in nuclear physics. The review by the Bochum 
group [13] summarizes the status of the field in 1978. 
From there on, applications in nuclear physics were more 
of a sporadic nature [2^. This was most probably due 
to the difficulty of hard NN interactions and their strong 
short-range repulsion and short-range tensor force. Mi- 
haila and Heisenberg employed coupled-cluster theory in 
the late 1990s [l3|- Their work culminated in the precise 
computation of the electron scattering form factor for 
^^O based on the Argonne vis potential combined with 
leading contributions from 3NF fld\ ■ 

Parallel to the field of nuclear physics, coupled-cluster 
theory saw its own career in ab-initio quantum chem- 
istry. After the pioneering works by Cizek psi, [s^l , the 
theory has become one of the main workhorses in quan- 
tum chemistry [13, [13, [H, [13] ■ The sheer number of 
applications and developments in that field de-facto es- 
tablished a "standard" or "canonical" way for how the 
method is being used to solve quantum many-body prob- 
lems. Reference [22| gives a summary of state-of-the-art 
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coupled-cluster calculations in quantum chemistry. 

Recently, coupled-cluster theory has seen a renais- 
sance in nuclear physics starting with the calculations 
of Ref. Tg"]. This approach differs from the one by Mi- 
haila and Hciscnberg as it employs coupled-cluster theory 
in the spirit of quantum chemistry and uses softer inter- 
actions. So far, the present approach has employed G 
matrices for the description of ground and excited states 
in ^He IITI and ^^O [l^, and for nuclei in the vicinity 
of ^^O [131 • The most recent calculations are based di- 
rectly on low- momentum interactions T^ow fc [H, , and 
the method has been developed to describe weakly bound 
and unbound helium isotopes within a Gamow-Hartree- 
Fock basis ^ and to include 3NF [2l|. 

Within coupled-cluster theory, the ground state of a 
mass A nucleus is written as 



(1) 
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and 
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Here and in the following, i, j, fc, . . . label occupied single- 
particle orbitals (as defined by the product state \(j))) 
while a,b,c, . . . refer to unoccupied orbitals. 

We take the reference state {(j)) as our vacuum state 
and normal-order the Hamiltonian H with respect to this 
state. In practice, we restrict ourselves to the truncation 
T = Ti + T2. This is the CCSD approximation, and the 
coupled-cluster equations are given by 
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is a np-nh exci- 
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is the similarity-transformed Hamiltonian (note that H 
is non-Hcrmitian) . The last expression on the right- 
hand side of Eq. ([7]) indicates that only fully connected 
diagrams contribute to the construction. The CCSD 
Eqs. ([5|) and ([6]) determine the amplitudes t° and tf^ 
of the Ip-lh and the 2p-2h excitation cluster operators, 
respectively. Once these nonlinear equations are solved, 
the amplitudes can be inserted into Eq. ([4]) to determine 
the ground-state energy. 

We remind the reader that an exact solution of the 
many-body problem would require us to employ the full 



excitation operator Eq. ([2]). Such a calculation is as ex- 
pensive as a full diagonalization of the Hamiltonian, and 
therefore impossible for medium-mass nuclei. CCSD is 
very efficient in the sense that it is a highly accurate ap- 
proximation with the investment of a modest numerical 
effort that scales as 0{n^n'^) with the number Uo of oc- 
cupied and the number Uu of unoccupied single-particle 
orbitals, respectively. The inclusion of the 3p-3/i clus- 
ter operator T3 would further increase the accuracy of 
the method. However, such CCSDT calculations come 
at the expense 0{n^n^) and, at present, are already pro- 
hibitively expensive compared to CCSD. For this reason, 
there is need for more approximate treatments of the full 
triples equations. 

There are various approximations to the full CCSDT 
equations, and the most popular of these schemes is the 
CCSD(T) approach [3l|. CCSD(T) includes diagrams at 
the CCSDT level that appear up to fifth order in per- 
turbation theory. It is a non-iterative approach since 
typically converged singles and doubles excitation ampli- 
tudes are used in the calculation of the triples energy 
correction. The CCSD(T) approximation is relatively 
inexpensive compared to CCSDT; no storage of triples 
amplitudes is required and the computational cost is a 
non-iterative 0{n^n'^) step. There is also a family of 
iterative triples correction schemes known as CCSDT- 
n 34]. Their derivation is based on perturbation theory 
arguments, 
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Here, F denotes the Fock operator (the one-body opera- 
tor that results from the normal ordering of the Hamilto- 
nian) . All these approaches require the storage of the full 
triples amplitudes t"^*^ and are therefore computation- 
ally considerably more expensive than the CCSD(T) ap- 
proach. However, for cases where the CCSD(T) scheme 
breaks down, one expects the CCSDT-n approaches to 
perform better. The latter approaches treat the triples 
corrections self-consistently and also involve the correc- 
tions 



(FT. 



(9) 
(10) 



to Eqs. (O and respectively. These corrections thus 
modify the values of the amplitudes tf and tf^. 



B. Low-momentum interactions and model spaces 

Nuclear interactions depend on the resolution scale at 
which details are probed and resolved. This resolution 
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scale dependence is similar to scale and scheme depen- 
dences in parton distribution functions. As a result, nu- 
clear interactions are defined by an effective theory for 
NN, 3N, and many-nucleon interactions and correspond- 
ing effective operators, 



F = 1/nn(A) + V^3n(A) + 



(11) 



where the momentum cutoff A denotes the resolution 
scale. Conventional nuclear forces are "hard" in the 
sense that they have large cutoffs that complicate few- 
and many-body calculations. These difficulties arise from 
high momenta and associated strong short-range repul- 
sion and short-range tensor forces, which lead to slow 
convergence with increasing basis size and requires re- 
summations in practice. 

Low-momentum interactions Mow k with variable mo- 
mentum cutoffs show great promise for nuclei (3^ . Issl . 
m, [si, [13, H H SI. Changing the cutoff leaves low- 
energy NN observables unchanged by construction, but 
shifts contributions between the potential and the sums 
over intermediate states in loop integrals. These shifts 
can weaken or largely eliminate sources of nonperturba- 
tive behavior such as strong short-range repulsion and 
the tensor force [H, [4lj . An additional advantage is that 
the corresponding 3N interactions become perturbative 
at lower cutoffs 3a] and are thus tractable in coupled- 
cluster theory f2ll| . The renormalization group (RG) evo- 
lution is implemented by coupled RG equations in mo- 
mentum space 14211 or by an equivalent Lee-Suzuki trans- 
formation [3, l45l |. 

The evolution to low- momentum interactions Mowfc 
weakens off-diagonal coupling and decouples the low- 
energy physics from high-momentum details [ii, [4^ . As 
a result, few- and many-body calculations converge more 
rapidly for lower cutoffs, which is important for extend- 
ing ab-initio approaches to heavier systems. Finally, the 
cutoff variation can provide estimates for theoretical un- 
certainties, which will be left to future work. In this 
paper, we use a sharp cutoff A = 1.9 fm^^ for the and 
"^He calculations, and Viow k is derived from the Argonne 
Wis potential [i^ in order to benchmar k ag ainst the Fad- 
deev and Faddeev-Yakubovsky results [33. For ^^O and 
""^Ca, we use a cutoff A = 2.1 fm^^ and compare to the 
importance-truncated NCSM study [23i] . 

Coupled-cluster theory is employed in a single-particle 
basis, and we use a model space consisting of spherical 
harmonic-oscillator states. The basis parameters are the 
number of orbitals and the oscillator frequency hu. Our 
largest model spaces include about 10^ single-particle or- 
bitals. 



III. RESULTS FOR AND *HE 

In this section, we present our coupled-cluster calcula- 
tions for the ground-state energies of and ''He, and we 
compare our results to the exact Faddeev and Faddeev- 
Yakubovsky energies of Ref. |35| . We first discuss in de- 
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FIG. 1: (Color online) CCSD and CCSD(T) energies for 
using a model space with fixed maximum N = 2n + I = 12 
and fixed hcj — 14 MeV as a function of the maximum orbital 
angular momentum I. For comparison, we also show the exact 
Faddeev result of Ref. . 
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FIG. 2: (Color online) CCSD and CCSD(T) energies for ''He 
as a function of the maximum orbital angular momentum I 
and the exact Faddeev-Yakubovsky (FY) result of Ref. [ssl ]. 
For details, see the caption to Fig. [T] 



tail the dependence on the size of the model space and 
the single-particle basis. The coupled-cluster calculations 
initially used a single-particle basis of oscillator states 
whose principal and angular momentum quantum num- 
bers n and / obey 2n -I- Z ^ A^, so -I- 1 is the number 
of major oscillator shells included. Note that in previ- 
ous calculations [3, [l^ [H, [l3| A^ denoted the number 
of major oscillator shells. However, we observed that the 
convergence with respect to the angular momentum I is 
much quicker, since only low partial waves contribute to 
low-energy properties, while the convergence with respect 
to the principal quantum number n is slower. This slower 
convergence is due to the sharp momentum cutoff used 
for Viowk- It is intuitively clear that a harmonic-oscillator 
representation of an interaction with a sharp cutoff needs 
a considerable number of radial wave functions to be ac- 
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FIG. 3: (Color online) CCSD and CCSD(T) results for the 
ground-state energy of as a function of the model-space 
size A'^ = 2n + I, with fixed I ^ 5 and fixed huj = 14MeV. 
For comparison, we also show the exact Faddeev result of 
Ref. m. 
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FIG. 4: (Color online) CCSD and CCSD(T) results for the 
ground-state energy of ''He as a function of the model-space 
size A'^ = 2n + I and the exact Faddeev- Yakubovsky (FY) 
result of Ref. [35l |. For details, see the caption to Fig. |3] 

curate. The recent work of Refs. [10, confirms this 
picture and demonstrates that smooth cutoffs improve 
the convergence in few-body calculations. 

Figures[T]and[2]show the convergence of our CCSD and 
CCSD(T) energies for ■^H and ''He using a model space 
with fixed N ^ 2n + I = 12 and fixed hj = 14MeV as 
a function of the maximum orbital angular momentum 
I. This implies that for Z = we include oscillator func- 
tions with n ^ 6 nodes; for I ^ I we include oscillator 
functions with n ^ 6 for the s states {I = 0) and n ^ 5 
for the p states {I = 1), and so on. Clearly, the angular 
momentum quantum number needs not to exceed / = 5 
for the ground-state energies of s-shell nuclei. Therefore, 
we limit our single-particle basis to I ^ 5 for the following 
coupled-cluster calculations of '^H and ^He. 

In Figs.[3]and[l we present our CCSD and CCSD(T) 
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FIG. 5: (Color online) CCSD and various approximate 
CCSDT energies for ^He as a function of the oscillator fre- 
quency fkj, for fixed A'' = 12 and 1^5. For comparison, we 
also show the exact Faddeev- Yakubovsky (FY) result. 

results for the ground-state energies of and ^He as a 
function of the model-space size N, with fixed I ^ 5 and 
fixed hw = 14 MeV. Both CCSD and CCSD(T) energies 
converge with respect to the model space size for N « 
12 . . . 14. For the largest model space with N = 16, we 
obtain for "^H, 

i^ccsofH) = -8.09 MeV, 
^^ccsD(T)fH) = -8.40 MeV, 

and for ''He, 

^^ccsdC^Hc) = -28.92 MeV, 
-£^CCSD(T) C'He) = -29.18 MeV. 

The CCSD(T) energies are within 70keV and lOkeV of 
the Faddeev and Faddeev- Yakubovsky (FY) results [sHi 
^(^H) ^ -8.470(2) MeV and £;('*He) = -29.19(5) MeV. 

Finally, wc study the dependence of our results on the 
oscillator frequency huj. This is shown in Fig. [5] for fixed 
TV = 12 and I ^ 5. While the CCSD results exhibit a 
very small variation over the shown Huj range, the vari- 
ation of the perturbative triples corrections CCSD(T) is 
somewhat larger. Moreover, the downward trend of the 
CCSD(T) energies with decreasing fku indicates that per- 
turbative triples corrections are starting to break down 
for smaller values of fko. The non- iterative perturbative 
triples correction assumes that we work in a basis where 
the Fock matrix is diagonal. However, our oscillator basis 
does not diagonalize the Fock matrix, so strict calcula- 
tions would have to iterate triples corrections until self- 
consistency is reached. From Fig. [5] we observe that iter- 
ative CCSD(T) improves on the non-iterative CCSD(T) 
results, but also has a downward trend with decreasing 
Huj. Finally, we present calculations based on the itera- 
tive CCSDT-1 approximation to fuU CCSDT. CCSDT-1 
includes all diagrams through fourth order in perturba- 
tion theory, but contrary to the perturbative CCSD(T) 
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FIG. 6: (Color online) CCSD and CCSD(T) results for the 
binding energy of using a model space with fixed maxi- 
mum iV = 2n + Z = 10 and fixed hio — 20 MeV as a function 
of the maximum orbital angular momentum I. 
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FIG. 7: (Color online) CCSD results for the binding energy of 
^® O as a function of the oscillator frequency hcu for increasing 
sizes of the model space N = 2n + I with fixed I ^ 7. 
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corrections, the CCSDT-1 approximation is treated self- 
consistently and the singles and doubles amplitudes are 
modified by the triples amplitude according to Eqs. ^ 
and We clearly find that CCSDT-1 improves on the 
triples corrections and leads to a very weak dependence 
on huj. Note that these CCSDT-1 results are also the first 
step towards full CCSDT calculations in nuclear physics. 

Our results for the light nuclei and *He demon- 
strate that coupled-cluster theory meets the benchmarks 
set by exact methods. It is therefore interesting to use 
this method to establish benchmark energies for heavier 
nuclei that other ab-initio approaches can compare to. 
This is the subject of the next section. 



IV. RESULTS FOR ^^^O AND *°CA 
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FIG. 8: (Color online) CCSD results (taken at the hco minima) 
for the binding energy of ^''O as a function of the model-space 
size N = 2n + I a.t fixed I ^ 7 and exponential fit (solid line) . 



Next, we present our coupled-cluster calculations for 
^^O and ^°Ca. In Fig.[6l we show the convergence of the 
CCSD and CCSD(T) results for the binding energy of 
"'^^O using a model space with fixed N — 2n + l — 10 and 
fixed hw = 20 MeV as a function of the maximum orbital 
angular momentum I. For ^^O, we find that I ^ 7 is suffi- 
cient to reach convergence at the 10 keV level. Therefore, 
we restrict the following coupled-cluster calculations for 
1^0 to Z < 7. 

Figure [7] shows the dependence of the CCSD binding 
energies of ^^O on the oscillator frequency hu for increas- 
ing sizes of the model space N with fixed I ^ 7. The 
largest calculations for = 13 include more than lO'^ 
single-particle orbitals. We observe that the CCSD ener- 
gies are converged at the 0.5 MeV level and can be used 
to extrapolate to infinite model space. This is demon- 
strated in Fig. [8] where we give the CCSD energies (taken 
at the huj minima) as a function of the model-space size 
N at fixed I ^ 7. Using the CCSD minima, we make an 
exponential fit of the form E{N) — Eoo + aexp{—bN) 



to the data points. The result is also shown in Fig. [5] 
and yields the extrapolated infinite model space value 
i^ccSD,oo(^^0) = -142.78 MeV. Our largest = 13 re- 
sult is 'i;ccSD(^^0) = -142.40 MeV. The conservative 
error estimate due to the finite size of the model space is 
thus about 0.5 MeV. 

In Fig. [21 we study triples corrections to the binding en- 
ergy of ^^O via the CCSD(T) and CCSDT-1 approaches 
as a function of the oscillator frequency hio for increasing 
sizes of the model space N {I ^7). We present results up 
to iV = 7 (eight major oscillator shells), which was the 
largest model space we could handle for the CCSDT-1 
scheme due to storage limitations. The CCSD(T) and 
CCSDT-1 energies agree nicely for the range of oscilla- 
tor frequencies and model spaces considered. The only 
difference is that, as for ^He, the CCSD(T) approach 
gives slightly more binding than CCSDT-1 for the largest 
model spaces. The close agreement between CCSD(T) 
and CCSDT-1 gives us confidence that the perturbative 
triples corrections work well over this regime. 
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FIG. 9: (Color online) CCSD(T) and CCSDT-1 results for the 
binding energy of as a function of the oscillator frequency 
hoj for increasing sizes of the model space N = 2n (I ^1). 
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FIG. 10: (Color online) Contributions to the binding en- 
ergy of ^''O from triples corrections CCSD(T) and CCSDT-1 
as a function of the model-space size N = 2n -\- I at fixed 
fujj = 22 MeV and I ^ 7. The dashed and dotted lines are ex- 
ponential fits and yield the extrapolated energy corrections. 
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FIG. 11; (Color online) CCSD results for the binding energy 
of ''"Ca as a function of the oscillator frequency hiu for in- 
creasing sizes of the model space N — 2n + I. (Note that 
there is no restriction on / for these model spaces.) 
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FIG. 12: (Color online) CCSD results (taken at the hw min- 
ima) for the binding energy of^OCa as a function of the model- 
space size A'^ = 2n+l (without restriction in I) and exponential 
fit (solid line). 



Let us study the contributions of the triples amphtudes 
to the binding energy of ^^O in more detail. Figure [TOl 
shows the energy differences AE — E — Eccsd that are 
due to triples corrections "(T)" and "T-1" as a function of 
the model-space size N at fixed huj — 22 MeV and / ^ 7. 
The corresponding exponential extrapolations to infinite 
model spaces yield —5.45 MeV from the "(T)" correction 
and —5.00 MeV from the "T-1" correction. This suggests 
that the error estimate for the CCSD(T) and CCSDT- 
1 calculations is about 0.5 MeV. This can be viewed as 
an error estimate due to the truncation of the cluster 
operator. Combined with the 0.5 MeV uncertainty due 
to the size of the model space, we thus arrive at an error 
estimate of about 1 MeV for ^^O. 

We note that "'^^O is overbound by about 20 MeV when 
compared to the experimental binding energy. A similar 
result has been found by Fujii et al. [i^ . However, a 
comparison of results based only on NN interactions to 



experiment is meaningless, since 3NF are crucial to de- 
scribe few- and many- body observables (see for instance 
the discussion in Refs. [2l|, [H, [11] ) . In nuclear matter, 
the corresponding 3NF contribution is repulsive and the 
expectation values remain consistent with chiral effective 
field theory power-counting estimates [36| . 

Finally, we turn to the more challenging case of '*°Ca. 
Figure [TT] shows the CCSD binding energy of ^"Ca as a 
function of the oscillator frequency hw for model spaces 
up to TV = 8 (nine major oscillator shells). (Note that 
there is no restriction on I for these model spaces.) This 
represents the largest coupled-cluster calculation to date 
in nuclear physics. In these largest calculations, we have 
40 active particles in 660 single-particle orbitals. The 
effective shell-model dimension in this space would be 
of the order of 10^^. From Fig. [Til we find that the 
CCSD energies of ''"Ca are converging reasonably well. 
We again expect that low-momentum interactions with 
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FIG. 13: (Color online) Contributions to the binding energy 
of ^''Ca from triples corrections CCSD(T) as a function of the 
model-space size N — 2n + I at fixed — 22 MeV (without 
restriction in I) and exponential fit (dashed line). 



smooth cutoffs will lead to even improved convergence. 
In Fig. [121 we present the CCSD energies for ^°Ca (taken 
at the fuv minima) as a function of the model-space size 
N. The exponential extrapolation to infinite model space 
yields -EccsD,oo('"'Ca) = -492.6 MeV, and we find that 
the CCSD energy for = 8 is about 4 MeV from the 
fully converged CCSD result. 

We then perform CCSD(T) calculations and show in 
Fig. [T3] the triples energy corrections as a function of the 
model-space size A'' at fixed Hlo = 22 MeV. Due to mem- 
ory limitations, we were not able to perform CCSDT- 
1 calculations in model spaces reaching up to iV = 7. 
The exponential extrapolation to infinite model space 
yields —11.70 MeV, while the largest N = 7 result is 
— 10.21 MeV. The convergence we find with respect to 
the size of the model space is similar to the recent results 
by Fujii et al. [50]. Recall that the different triples cor- 
rections for ^^O differed by about 10% from each other. 
Thus, we estimate that the error due to the truncation 
of the cluster operator is about 1 MeV for "^^Ca. The 
total error estimate for "^^Ca is thus about 5 MeV, and 
is dominated by the uncertainty due to the finite size of 
the model space. 

We summarize our coupled-cluster results for the bind- 
ing energies of ^He, ^^O, and ""^Ca in Table |T1 which 
gives the extrapolated correlation energies A£^ccsd and 
A£;ccSD(T)- We find that for ^He, ^^O, and ^'^Ca, the 
triples corrections are a factor of « 0.015,0.066, and 
0.081 smaller than the CCSD correlation energies. From 
this, we again estimate the missing correlation energy 
from quadruples, pentuplets, and so on, to be of the or- 
der of 1 MeV for ''"Ca. We note that ^^O is overbound 
by about 20 MeV and ^°Ca by about 150 MeV when com- 
pared to the experimental binding energies. This is not 
surprising and points to the importance of 3NF for nu- 
clear structure calculations [2l|, IH, . 

We can compare the coupled-cluster energies to the 
recent importance-truncated NCSM results of Roth and 





^He 


16q 


«Ca 


Eo 


-11.8 


-60.2 


-347.5 


A-BccsD 


-17.1 


-82.6 


-143.7 


^-£'CCSD{T) 


-0.3 


-5.4 


-11.7 


-E'CCSD{T) 


-29.2 


-148.2 


-502.9 


exact (FY) 


-29.19(5) 







TABLE I: Reference vacuum energies, Eo, CCSD and 
CCSD(T) correlation energies, A_Eccsd and Ai?ccsD(T), and 
binding energies -Eccsd{t) for *IIe, ^^O and ^'^Ca. The vac- 
uum energies, -Eo, are for hu) — 14 MeV in the case of ''He and 
huj = 22 MeV for ^''O and *°Ca. The CCSD and CCSD(T) 
energies are the extrapolated infinite model space results. The 
exact Faddeev-Yakubovsky result is from Ref. [35i ] 



Navratil [23| which are based on the same Mowfc inter- 
action. The importance-truncated NCSM combines a 
particle- hole truncation scheme (4p-4/i for ^^O and 3p-3h 
for ^"Ca in Ref. [23]) with importance sampling of many- 
body states based on perturbation theory. The particle- 
hole truncation scheme leads to unlinked diagrams and 
hence is not size-extensive [2^ . Using an exponential 
extrapolation. Roth and Navratil [1^ find binding ener- 
gies E = -137.75 MeV for ^^O and E ^ -461.83 MeV 
for '^''Ca at the minimum in tiuj. Our coupled-cluster 
results idicate that the converged energies are approxi- 
mately 10 MeV and 40 MeV lower for and '^'^Ca, re- 
spectively. Note also that CCSD scales computationally 
more favorably than a full Ap-Ah calculation, while it al- 
ready includes a considerable part of linked Ap-Ah exci- 
tations il. 



SUMMARY 



In summary, we have performed ab-initio coupled- 
cluster calculations for ^H, "^He, ^^O, and ''"Ca based 
on low- momentum interactions Mowfc- At the CCSD(T) 
level, the ground-state energies for and ''He are prac- 
tically converged with respect to the size of the model- 
space and exhibit a very weak dependence on the oscil- 
lator frequency. The resulting energies are within 70 keV 
and lOkeV of the corresponding Faddeev and Faddeev- 
Yakubovsky benchmarks. For ^^O and ''"Ca, we estimate 
that the CCSD(T) binding energies are converged within 
1 MeV and 5 MeV, respectively. Future calculations will 
include convergence studies for low-momentum interac- 
tions with smooth cutoffs [H, |4^ and advancing the 3NF 
frontier to medium-mass nuclei based on the findings of 
Ref. [2l|. Our results confirm that coupled-cluster the- 
ory is a powerful ab-initio method that meets and sets 
benchmarks. 
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